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ABSTRACT 

We use a simple model of spheroid formation to explore the relationship between the cre- 
ation of stars and dust in a massive proto galaxy and the growth of its central black hole. This 
model predicts that sub-mm luminosity peaks after only ~ 0.2 Gyr. However, without a very 
massive seed black hole, Eddington-limited growth means that a black hole mass of 10 9 M , 
and hence very luminous AGN activity, cannot be produced until > 0.5 Gyr after the forma- 
tion of the first massive stars in the halo. The model thus predicts a time-lag between the peak 
of sub-mm luminosity and AGN luminosity in a massive proto-elliptical of a few times 10 s 
years. For a formation redshift z ~ 5, this means that powerful AGN activity is delayed until 
z ~ 3.5, by which time star formation in the host is ~ 90% complete, and sub-mm luminosity 
has declined to ~ 25% of its peak value. This provides a natural explanation for why success- 
ful sub-mm detections of luminous radio galaxies are largely confined to z > 2.5. Conversely 
the model also predicts that while all high-redshift luminous sub-mm-selected sources should 
contain an active (and growing) black hole, the typical luminosity of the AGN in such objects 
is ~ 1000 times smaller than that of the most powerful AGN. This is consistent with the al- 
most complete failure to detect sub-mm selected galaxies with existing X-ray surveys. Finally 
the model yields a black-hole: spheroid mass ratio which evolves rapidly in the first Gyr, but 
asymptotes to ~ 0.001 — 0.003 in agreement with results at low redshift. This ratio arises not 
because the AGN terminates star formation, but because fueling of the massive black hole is 
linked to the total mass of gas available for star formation in the host. 

Key words: galaxies: elliptical and lenticular, cD — galaxies: active — galaxies: evolution 
— quasars: general — stars: formation — (ISM:) dust, extinction — submillimetre 



1 INTRODUCTION 

Until recently the formation/evolution of active galactic nuclei and 
the formation/evolution of galaxies of stars have been viewed and 
studied as separate problems in observational cosmology. However, 
two recent observational discoveries have changed this. First, it is 
now clear that the global history of star-formation and AGN ac- 
tivitv_track£arl^otiiej^athe^ 



lop |!997| ; |Boyle & Terlevich||l998| : jpranceschini et al|[200C| ). Sec- 
ond, it is also now well established that essentially all present- 
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day galaxies contain a central black hole whose mass is propor- 
tional to the mass of the sph eroidal component of the ho st, with 
M bh ~ 0.0 01-0.003 M sr>h (iKormcndy & Gcbhardt||200l|; Merritt 
& Ferrarese 2001). Consequently the growth and fueling of super- 



massive black holes and the construction of galactic spheroids are 
now se^nasjjotovh^dlyinrim & 



Rees |1998| ; ^abian||l999| ; pranato et alj|200l| ). 

This increased awareness of possible links between black-hole 
and spheroid formation has coincided with recent observational ad- 
vances which have the potential to provide greatly-improved con- 
straints on the evolution of massive galaxies and AGN. In particular 
new sub-mm surveys with the SCUBA camera on the JCMT have 
resolved a substantial fraction of the sub-mm background (Smail 
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et al. |1997| |Hugheset al||l998| ; pales et alj|l999| ), and have revealed 
a population of massive, dusty star-forming objects at z > 1.5 — 2 
which are sufficiently numerou s to account for the formation o f all 
present-day massive ellipticals ( Dunlop] 200 1 ; 



Scott et al. 



2002). At 



the same time, deep X-ray surveys have resolved the X-ray back- 
ground and have enabled the evo lution of X-ray em itting black 
holes to be traced out to z ~ 4 (Miyaji et al. 2000). Moreover, 



motivated by the potential links between black-hole and spheroid 
formation outlined above, there is considerable current interest in 
determining the level of overlap between the sub-mm and X-ray 
populations. Current obser vations indicate this overlap is, perhaps 
surprisingly, small 
Barger et al][200l|) 



10% (Fabian et al. 



200C ; Almaini et al 



2002 



An alternative route to exploring the link between black-hole 
and spheroid formation is provided by observations designed to 
measure the star-formation history o f the hosts of known a ctive 
black holes as a function of redshift. Archibald et al. (2001) per- 
formed this experiment for the most massive black holes by un- 
dertaking a sub-mm study of a sample of ~ 50 radio galaxies 
spanning a wide range in redshift 1 < z < 5 at constant radio 
power. At low redshift it is now well established that powerful ra- 
dio sources are produced by only the most massive black holes 
with Mbh > 10 9 Mpj, residing in mas sive elliptical hosts with 
M aph ~ 1O 12 M ( fc)unlop et at|[2002| ). Thus, sub-mm observa- 
tions of powerful radio galaxies over a wide range of redshifts can 
offer insight into the formation history of at least a subset of the 
most massive elliptical galaxies. 

The principle result of this study was that the sub-mm de- 
tectability of powerful radio galaxies was found to be a strong func- 
tion of redshift, with a detection r ate of ~ 75% at z > 2.5 com- 
pared to only ~ 15% at z < 2.5 (Archibald et al. 2001). Indeed, 



the average sub-mm luminosity of the radio galaxies in the sam- 
ple could be approximated as rising with increasing redshift with 
L oc (1 + z) 3 out to z > 3. 

By definition, all the objects observed in this study contained 
sufficient gas to fuel the activity of their central black holes at the 
epoch of emission. Viewed in the light of this selection bias, the ob- 
served strong redshift dependence of sub-mm luminosity provides 
rather convincing evidence that the hosts of massive black holes 
at high redshift are very different to the relatively passive elliptical 
hosts of low-redshift AGN. Specifically, irrespective of whether the 
dust is heated primarily by the UV output of young stars or by the 
AGN itself, these results indicate that the mass of dust (and hence 
gas) in high-redshift AGN hosts is much greater at z > 2.5. The 
obvious implication is that by z ~ 3 we are probing an era in which 
a significant fraction of the eventual stellar population of a typical 
massive elliptical has yet to be formed. 

Interestingly, a comparison of the redshift distribution of the 
sub-mm detected radio galaxies with current best-estimates for the 
redshift distribution of sub-mm s elected sources indica t es that they 



are st atistically indistinguishable ( Archibald et al. |200l| 



Smail et al 



200C). This provides some justification for believing that the sub- 
mm evolution observed for the radio-galaxy sample may apply to 
spheroid formation in general. However, this then raises the obvi- 
ous counter question of why so few bright sources discovered via 
sub-mm surveys display any obvious signs of powerful AGN activ- 
ity? 

Motivated by these new observational results we here develop 
further a simple model for the dust-enshrouded formation of mas- 
sive spheroids which we have previously applied to explore the 
red envelope of galax y evolution at optical-infrared wavelengths 
(limenez et al. 1999). This model is constructed by combining 



t he ste l lar po pulation and dust-formation models of Jimenez et al 
(2000, 2002) with a multi-zone chem o-hydrodynamic model fo r 
a collapsing gas spheroid developed by Friaca & Terlevich (1998). 



The model yields predictions for the time evolution of gas, dust and 
the stellar population of the galaxy as a whole, but also allows us to 
set constraints on the growth and fueling of a central black hole by 
tracking the rate at which gas is deposited into the central 100 pc of 
the galaxy. The predictions of this model not only provide insight 
into the results of the SCUBA radio-galaxy survey, but also offer 
a possible explanation for the relative failure of X-ray surveys to 
detect bright SCUBA sources. The model also demonstrates that a 
correlation between black-hole and spheroid mass consistent with 
that seen in low-redshift galaxies can be produced without invok- 
ing a direct causal connection between black-hole emission and the 
term ination of star formation in the host (as suggested by S ilk & 



Rees |1 998| and |Fabian| |l 99sj )~ 



At first sight, it might appear that this spherically symmetric 
calculation of spheroid formation within an isolated and complete 
dark matter halo is of dubious relevance in the context of the ac- 
cepted paradigm of hierarchical galaxy formation driven by halo 
mergers. However, to date, multi-zone chemo-hydrodynamic mod- 
els of the type considered here have only been developed in one 
dimension. Despite this technical limitation, we feel it is neverthe- 
less important to consider the detailed predictions of such a model, 
especially given the extremely simplified recipes of star-formation 
currently adopted in t he construction of semi-anal y tic models of 
galax y formation(e.g, |Somerville & Primack 



1999; Benson et al 



2001). Moreover, for AGN host galaxies forming at z ~ 3 (i.e. at 



or near the peak epoch of quasar activity), the simple model pre- 
sented here may in fact offer a good approximation to the average 
behaviour of the baryons within the forming host galaxy. This is be- 
cause, in such highly-biased regions within the early Universe, the 
merging of collapsed dark matter haloes of mass ~ 5 x 10 10 M 
into a final halo mass ~ 10 13 Mq is completed very rapidly (see 
Section ||). 

The structure of this paper is as follows. In Section |^ we de- 
scribe the key components of our theoretical model of a young 
elliptical galaxy, and summarize the predicted time evolution of 
the observable properties of such an object. The predictions of this 
model at sub-mm wavelengths are then compared with the results 
of the sub-mm study of radio galaxies in Section^. In Section ^ we 
explore the predicted fueling and growth of the central black hole, 
and compare the predictions of the model with measurements of 
black-hole masses at low redshift and with the cosmological evolu- 
tion of AGN activity. Finally, in Section ^ we discuss the implica- 
tions of this model for the connection between black-hole and host- 
galaxy evolution, focusing on the predicted sub-mm detectability of 
AGN, and the X-ray detectability of sub-mm selected sources. For 
the calculation of physical quantities we assume a flat universe with 
Q.M = 0.3, Q.a = 0.7 and H = 70 kms^Mpc" 1 throughout. 



2 MODEL OF A FORMING ELLIPTICAL 

The stellar populations observed in massive elliptical galaxies, es- 
pecially those residing in clusters, appear to be highly coeval, con- 
sistent with the majority of the star s having formed on a timescale 

Whether these stars formed 



1 Gyr (e.g. 



Bower et al 



in a single region, or in smaller separate regions which merged to 
form the final spheroid, the star formation is expected to be vigor- 
ous, comprehensive, and obscured by dust - the initial gas reservoir 
being quickly consumed in the star formation process, and the rem- 
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nants of the process being blo wn away by stellar winds and super- 



novae (e.g., Zejrf&SnkJ 



1996) 



Jimenez et al. ( 200C ) have developed multi-frequency models 
of galaxies formed via such rapid and violent dust-enshrouded star- 
bursts. These models predict the time evolution of the ultraviolet- 
millimeter spectrum of a young elliptical galaxy based on a calcula- 
tion of the star-formation rate and the resultant chemical evolution 
of the interstellar medium (ISM). 



Event 


Timescale 


50% stars formed 


0.3 Gyr 


98% stars formed 


1.0 Gyr 


100% stars formed 


1.7 Gyr 


dust mass peaks 


0.2 Gyr 


metallicity peaks 


2.0 Gyr 



Table 1. Key stages of formation as predicted by the dusty starburst model. 



2.1 How the Models Work 

There are three key components to the model: (i) the stellar emis- 
sion from the model galaxy, (ii) the re-processing of this emission 
by dust, and (iii) the star-formation law which underpins both these 
processes. 



2.1.1 Evolution of the Stellar Population 

Given a star-formation law and an initial mass function, the stellar 
populations that for m in each time step build up to form a com- 
plete model galaxy (fimenez et al. 1998, 2002). Given a value for 



metallicity and the initial stellar mass (governed by the IMF), the 
evolutionary track of a star is calculated by the c ode JMSTAR15 
( Jimenez et al. 1996; limenez & MacDonalc 1996), which includes 



a proper treatment of the post-main-sequences phases of evolution. 
In the absence of dust and gas the integrated emission of the model 
galaxy at any time step is simply the sum of the flux from every star 
in the population. 

The model assumes a Salpeter IMF of (f>{M) oc M~<- 1+x \ 
where <f>(M) is the number of stars formed per unit volume with 
a mass M, and x = 1.35. The IMF includes a low-mass cutoff at 
0.1 Mq, and is assumed to be constant in both time and space. 



2.1.2 Dust Model 



the formalism of Xu & de Zotti 



Th e dust model ( limenez et al. 200C ) uses a simplified version of 
the Draine & Lee J 1984 ) dust emissio n template, and is based on 



(1989) 



For the purposes of calculating the far-infrared emission, the 
dust is simply assumed to lie in a thin shell which surrounds the 
entire galaxy, illuminated by the uniform radiation field generated 
by the stellar population. This distribution is not realistic. However, 
individual stars and gas clouds cannot be resolved in high-redshift 
systems. Given the lack of knowledge concerning the true distribu- 
tions, this simple model is no worse than a more complicated model 
of dust and star geometry which cannot possibly be justified by the 
data. 

This simple dust distribution is not, however, deemed ap- 
propriate for calculating the reddening of the galaxy light. Dense 
molecular clouds are the birth place of stellar clu sters, where all the 
massive stars are formed (e.g., |^ada et aL 1991). It takes 1.5 x 10 7 
years for stellar winds and supernovae to disperse these clouds, 
which is approximately the lifetime of the most massive stars (e.g., 

|Blitz & Shij |198C| ). The 



Herbig 



1962 



Elmegreen & Lada 



1977 



primary source of extinction of the galaxy light is the dust inside 
these molecular clouds, not the dust lying outside the starburst re- 
gion (which is the assumed geometry for estimating the far-infrared 
emission of the galaxy). Thus, the model assumes that for the first 
1.5 x 10 7 years of a star's life, the star is only visible in the infrared, 
with the optical and ultraviolet emission completely obscured. Af- 
ter 1.5 x 10 7 years, the molecular cloud is assumed to have been 



dispersed, and the star becomes visible in the optical and ultraviolet 
regimes. 

2.1.3 Star- Formation and Chemical Evolution 



The modjeluse^ji^tar-fonr^^ & Ter- 

levich ( 1998| ). Assuming spherical symmetry, the galaxy is divided 
into 100 radial zones. A 1-D hydro-dynamical code is used to fol- 
low the radial movement of the ISM and stars, including episodes 
of gas infall and outflow. The star-formation rate has a power- 
law dependence on the gas density: SFR(r,t) oc p 0,5 . The star- 
formation rate also has a radial dependence on the flow of the gas. 

Using the resulting distribution of the gas flow, the chem- 
ical evolution is solved by taking into account the lifetimes of 
stars, and using modern nucleosynthesis predictions to compute 
the yields of stars of different masses. In particular, for stars with 
masses ^Oijj^^Jl^A^^^^J^hcjidds^ 

& Vo li dl98l|), while fo r stars with 8 < M/M G) < 10 the yields 
from 
10 < 



Hillebrandt (1982) were assumed . For SNII progenito r s with 



M/Mq < 40 we adopted the Woosley & Weave: (1995) 



yields and extrapolated these for M > 40M/M a Finally, for SN 
la we adopted the yields from Nomoto et al ( 1984). 



As star formation continues the interstellar gas is heated by 
Type II supernovae. Eventually, the thermal energy of the gas is 
enough to overcome the escape velocity, and the remaining gas is 
expelled from the galaxy. At this point, the primary burst of star 
formation ceases. Later on, the red giants of low mass stars and 
Type I supernovae replenish the gas reservoir and metal reserves 
of the ISM, albeit to a limited extent. There is also a fresh infall 
of primordial gas which further dilutes the metallicity of the ISM. 
These new gas reserves allow residual star formation to occur. 



2.2 Time Evolution of the Models 

Friaca & Terlevich performed their 1 -dimensional hydrodynami- 
cal/chemico calculations for a range of assumed proto-galaxy gas 
masses, a fact which we exploit in Section 4.2. However, for com- 
parison with the properties of radio galaxies we adopt a fiducial 
model which yields an eventual stellar mass of ~ 10 12 M . The 
key stages of formation of a massive elliptical galaxy, as predicted 
by this dust-enshrouded starburst model are summarized in Table[I[ 

The primary burst of formation is rigorous and intense, tak- 
ing only ~ 0.3 Gyr to form the first half of the stellar popula- 
tion. At its peak, the star-formation rate is no higher than a few 
1000 Meyr -1 . As the primary burst of formation draws to a close, 
subsequent episodes of formation occur due to the replenishment 
of the gas reservoir. This residual formation is less vigorous, and 
although 98% of the stellar population is formed within a Gyr, it 
takes a further 0.7 Gyr to form the final 2%. 

One prediction of this model which appears to be a generic 
prediction of models of dust-enshrouded star formation is that the 
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dust mass peaks when the gal a xy is appro ximately half formed. Fo r 
example, Edmunds & Eales (1998) and Frayer & Brown (1997) 



have developed analytical equations which govern the relationship 
between gas depletion, star-formation, and dust creation. These 
equations also predict a maximum in dust mass half-way through 
the star formation process. This is a consequence of competing pro- 
cesses: star formation is required to create the dust, but star forma- 
tion also consumes the metal-enriched gas which houses the dust 
produced by previous generations of massive stars. It transpires 
that these competing effects balance to produce a maximum in dust 
mass approximately half-way through formation of the stellar pop- 
ulation. 

Even though half of the stars are formed very rapidly, in the 
first 0.3 Gyr, it can be seen from Table 1 that star formation persists 
in the host at some le vel for more than a Gyr. As demonstrated by 
Jimenez et al. (1999), at optical-infrared wavelengths this has the 



important effect of predicting that, even for a formation redshift 
Zf — 5, the largest values of R — K observed from an evolving 
elliptical are produced at z ~ 1.5 — 2. This prediction counters 
previous suggestions that pseudo-monolithic models for elliptical 
formation can be excluded in favor of hierarchical m odels simpl y 
due to lack of extremely red objects at high redshift (Zepi 1997). 



In the current context, this extended tail of star-formation activity 
leads to the prediction that a massive elliptical should be detectable 
at the 1-2 mjy level at 850^tm up to ~ 2.5 Gyr after the onset of 
star-formation activity. 



a. 
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3 GALAXY EVOLUTION - COMPARISON WITH 
RADIO-GALAXY DATA 

3.1 An Example: 4C41.17 

The multi-frequency evolution of galaxy luminosity predicted by 
the model can be used to assess the evolutionary status of well- 
observed massive galaxies at high redshift. As an example of how 
this can work we here consider two specific predictions of the 
model. The first is the predicted evolution of sub-mm flux density, 
and the second is the predicted evolution of optical-sub-mm color. 
The first of these is obviously an absolute property, dependent not 
only on the viability of the model but also on galaxy mass and as- 
sumed cosmology. In contrast the second is a relative measure of 
the mass of dust compared to the mass of stars which have already 
formed. While reliable estimates of galaxy mass at high redshift are 
not yet achievable, one can at least ask the question of whether, for 
reasonable values of assumed mass, these two predictions lead to 
the same conclusion about the evolutionary status of a given object. 

The z — 3.8 radio galaxy 4C 41.17 is one of the bes t studied 
of al l objects known at z > 3 (e.g., [chambers et al.||l99C|; Dunlop 



et al. |1994| ; |Dey et alj|1997| ; |van Breugel et al.||1998[ ). Two other 
properties of 4C41.17 make it particularly useful as a test of the 
model predictions. First, as the host of one of the most powerful 
known radio sources, it would seem safe to assume that it is the 
progenitor of a massive elliptical with an eventual stellar mass of at 
least 10 12 Mq . Second, its observed optical emission is known to 



be dominated by starlight rather than the AGN ( Dey et al 
In Figure [j], 

-£'850Mm/-I'2.2/jm/(l + z) 



1997) 



we compare 



and M, 



d lis t 



the time evolution of 
predicted by the model 
with the observed values for 4C41.17. This comparison indicates 
that 4C41.17 is experiencing a starburst which is ~ 0.5 — 0.8 Gyr 
old, having begun at z/ ~ 5 — 6. If the model is a correct 
description of 4C41.17, then ~ 75% of the final stellar mass of the 
galaxy has already been formed. 
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Figure 1. Evolution of Lsso^m/L 



pm/-tJ2.2/jm/(l+z) 



and M,i 



with time as 



predicted by the dusty starburst model (dashed lines). The corresponding 



bounds indicated by the shaded region ( 


Archibald et al. 


2001; 


etal. 


19%; 


<\rchibalc| 


19991. 



van Breugel 



This conclusion is driven primarily by the value of the ratio 



/L 2 .2|itn/(l+z) 



observed for 4C41.17 which, in essence, 
indicates that while active star-formation is still continuing at a high 
level ~ 1000 MQyr _1 the majority of the stellar population is al- 
ready in place. But we can see that this conclusion is also consistent 
with the observed 850/im flux-density of 4C41.17 provided its to- 
tal baryonic mass exceeds ~ 3 x 10 M©, which seems perfectly 
reasonable g iven the stellar mas ses of the hosts of low-redshift ra- 
dio sources (t>unlop et ai]|2002|). Thus, consistent with the original 
conclusions of Dunlop et al (1994), this model comparison indi- 
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Figure 2. Redshift evolution of a massive spheroid as predicted by the 
model for an assumed formation redshift zj = 5.0. The dotted line rep- 
resents the fractional stellar mass, the dashed line represents the fractional 
gas mass, and the solid line shows the dust mass normalized to its max- 
imum value. The data points give the weighted-mean 850- um luminosity 
<* f/cc n ,, .,, >. of the radio galaxy sample binned in redshift ( Archibald et al 
2001 ). The solid circles indicate the entire dataset; the opei 
a sub set chosen to remove any radio-luminosity bias (see |Archibald et al 
2001). 



cates that 4C41.17 is observed in the final stages of its formative 
starburst. 

If this conclusion is correct, it is worth noting that the model 
predicts an object as massive as 4C41.17 would have produced a 
sub-mm flux density as large as S&50 — 15 — 20 mjy if it had been 
observed only 0.3-0.6 Gyr earlier, only 0.2 Gyr after formation (i.e. 
at z = 4.5 for a formation redshift 2/ = 5). However, if observed 
at this stage in its evolution the model also predicts that the galaxy 
would have been much harder to detect in the near-infrared, with 
K ~ 21. Moreover, as discussed further below, at this stage it 
seems unlikely that 4C41.17 would have been capable of power- 
ful radio emission. The inferred properties of 4C41.17 only ~ 0.5 
Gyr prior to the actual observed emission epoch are therefore inter- 
estingly similar to the observed properties of many of the sub-mm 
select ed sources found in the blank-field SCUBA surveys (Dunlop 
2001 ). This point is pursued further in section 5. 



3.2 Redshift dependence of sub-mm luminosity 

Assuming a formation redshift of zf = 5.0 as inferred above for 
4C41.17, we plot in Figure the redshift evolution of the gas mass, 
dust mass, and stellar mass predicted by the model. 

The formation of the spheroid lasts approximately 1 Gyr and 
is ~ 95% complete by z ~ 2.5. While it would seem unreason- 
able to ascribe a single formation redshift 2/ = 5 to all power- 
ful radio galaxies, it is nevertheless interesting that 2 ~ 2.5 is 
also the redshift below which radio galaxies generally appear as 
relaxed spheroids in infrared images, and above which they gen- 
erally appear to be more complex, extended and clumpy sources 



in the rest-frame optical ( van Breugel et al 



1998 



Pentericci et al 



2001]). Moreover, as also illustrated in Figure [2| the predicted red- 
shift evolution of sub-mm luminosity is undeniably reminiscent of 
the redshift depend ence of mean Lgso^m for the full radio-galaxy 
sample observed by 



Archibald et al 



(2001) 
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In summary, the observed sub-mm - optical properties of high- 
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Figure 3. Top panel: predicted evolution of apparent /C-band magnitude as 
a function of redshift. Bottom panel: predicted evolution of absolute i?-band 
magnitude as a function of redshift. The dotted line shows the correspond- 
ing absolute R magnitude for a present-day L* galaxy. 



redshift radio galaxies are, at least to first order, consistent with the 
model predictions for an assumed typical formation redshift Zf — 
5. 



3.3 Redshift dependence of optical luminosity 

Finally, we plot in Figure |^ the redshift evolution of optical lumi- 
nosity, both in terms of absolute i?-band magnitude, and in terms 
of observed K magnitude. As can be seen from the lower plot, the 
2 = i?-band absolute magnitude is 1 .4 magnitudes brighter than 
M^, consistent with observations which indicate that radio galax- 
ies are typically 3 — 4 times more luminous than L* galaxies. 



4 BLACK HOLE EVOLUTION 

Having established that the model provides a plausible description 
of the evolution of massive ellipticals as traced by radio galaxies, 
we now explore the predictions of the model for central black-hole 
growth, and AGN output as a function of time/redshift. 
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4.1 Black Hole Growth 



10 F 



As detailed in Friaca & Terlevicl (1998) one output of the multi 



zone model is the amount of gas delivered per unit time to the cen- 
tral regions of the galaxy (r < 100 pc). 

To calculate the growth of the central black hole we assume 
that all of this fuel is available for consumption by the central black 



hol e (true in the absence of any other limiting factor, B urkert & 
Silk [2001 ). However, as illustrated in Figure |4| (for our fiducial 



model M ga i = 10 Mq) in practice for most of the first Gyr, the 
growth of the central black hole is limited not by this fuel supply, 
but by the Ed^mgtmiJmviKfoll^ & 
Blandford |1992[ ), unless a very large seed mass (~ 100000 Mq) is 
adopted. Figure ^a shows how black-hole mass consumption com- 
pares with central fuel delivery as a function of time for an assumed 
seed mass of 10 Mq and an assumed black-hole radiative efficiency 
of 10%, while Figure^] shows the effect of increasing the seed mass 
to 1000 Mq. We consider 10 — 1000 Mq to represent a reasonable 
range for assumed seed mass because obser yationally no stars with 
masses larger than ~ 200 Mq are known ( Figer et al. 1998 ), and 
because theoretical models indicate that the first generation of pop- 
ulation III stars may have produced blac k holes of mass ~ 100 Mq 



( Fryer et al 



2001 



Madau & Rees 2001) 



Figures |pandl|b show how black-hole mass grows as a func 



tion of time for both assumed seed masses. In the first case the final 
^ Mq , while in the second case it is 

-|9 i 



black-hole mass is 1.3 x 10 

1.6 x 10 9 Mq. Thus, varying the seed mass by a factor of 100 has 
only a very minor effect on the final mass of the black hole. Indeed, 
even if it is assumed that all the mass deposited into the central 
100 pc is ultimately consumed by the black hole (i.e. including that 
fraction of gas which is delivered, but cannot be consumed during 
the Eddington-limited accretion phase) the final mass is still not 
very different, at 2.0 x 10 9 Mq. 

Inspection of Figures ^ and ^| shows that the black hole grows 
and radiates at the Eddington limit for the first ~ 0.5 — 0.7 Gyr 
(doubling in mass every 3 x 10 r yrs) but that beyond that point, 
when the black hole has grown to a mass of a few times 10 8 Mq, 
the fuel delivery is unable to sustain further Eddington-limited fu- 
eling/growth. However, the black hole continues to grow to its final 
mass in excess of 10 9 M@ for a further Gyr or so, through sub- 
Eddington accretion before the gas supply is essentially switched 
off 3 Gyr after formation. 

Although this paper does not pretend to place any new con- 
straints on our understanding of black-hole accretion mechanisms, 
it is nevertheless important to ask whether the accretion rates pro- 
duced by this model are consistent with current understanding of 
the rate at which a black hole can consume matter from an ac- 
cretion disc. In particular, one might question whether the black 
hole is capable of consuming the gas delivered to it at the peak rate 
which approaches 2 Mgyr -1 (see Figure^). In fact, for a massive 
spheriod (with central velocity dispers ion a sp h — 200 kms -1 ) the 
calculations of Burkert & Silk (2001) indicate that accretion disc 



viscosity can be expected to limit the mass consumption rate of 
a super-massive black hole at the centre of a forming spheriod to 
~ 2 — 20 Mg yr -1 for a critical R eynolds number R cr in the range 
1000 — 100 ( Duschl et al. 200C). Thus, our current knowledge of 
accretion disc physics does not appear to conflict with our assump- 
tion that the black-hole can consume all the gas delivered to the 
central region of the galaxy at a peak rate which, in this particular 
model, approaches 2 Mgyr -1 . 

Several features of these results deserve comment. First, due 
to the assumption of Eddington-limited growth, the spheroid does 




Figure 4. Black-hole fueling rate as a function of time, assuming initial 
black-hole seed masses of 10 Mq and 1000 Mq . The solid line depicts the 
rate at which gas is deposited into the central 100 pc of the model galaxy. 
The dashed line represents the rate of gas consumption by the black hole 
assuming that this rate cannot be greater than that dictated by the Eddington 
limit for an assumed radiation efficiency of 10%. The shaded area indicates 
that Eddington-limited accretion persists until ~ 0.5 — 0.7Gyr, beyond 
which the central cooling flow cannot sustain Eddington-limited growth. 
Most of the final black-hole mass is then accreted at sub-Eddington rates 
over the subsequent ~ 1 Gyr. 



not contain a black hole of sufficient mass to host a quasar (i.e. 
M bh > 10 s M ; Dunlop et al. 2002) until 0.5 - 0.7 Gyr after 
commencement of star-formation. This predicted delay is still con- 
sistent with the discovery of quasars at z ~ 6 in our favored cos- 
mology (t un i Ver3e = 0.9 Gyr at z = 6), and can only be signif- 
icantly accelerated by the assumption of either a very large seed 
mass or possibl y runa way merging of black holes as suggested by 



Ebisuzaki et al 



(2001) 



Second, the final black-hole mass results in a black- 
holeispheroid mass ratio of 0.0011 - 0.0014 (or 0.0017 if all avail- 
able gas is consumed by the black hole), i n very good agreement 
with observational results at low redshift (Kormendy & Gebhardt 



2001 



Merritt & Ferrarese 2001 ). Note that this ratio arises not from 



invoking the assumption that the black-hole itself is responsible 
for r egulating and terminating star-formation in the host ( Silk & 



Rees |1998[ ^abian||l999| ), but ra her because the gas dynamics and 
star-formation in the forming spheroid has not allowed Eddington- 
limited exponential black-hole growth to continue beyond the first 
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Figure 5. As Figure Q, but this time plotting the the logarithm of black-hole 
mass as a function of time, for both alternative assumed seed masses. 



Gyr, and ultimately terminates significant black-hole fueling after 
the first 2 — 3 Gyrs (we note that Burkert & Silk 2001 also appeal 
to star formation to limit excessive black-hole growth, but in their 
model the star-formation in question is specifically assumed to be 
in the outer accretion disc). 

Third, there is a window lasting ~ 1 Gyr where the black hole 
has reached a mass sufficient to be capable of producing a powerful 
AGN and during which it is still receiving fuel at a sufficient rate 
(~ 1 Mgyr" 1 ) to do so. The bulk of the mass of the black-hole is 
actually accreted at sub-Eddington rates during this 'AGN epoch'. 

Fourth, the epoch of peak sub-mm luminosity precedes peak 
AGN luminosity by 0.3 — 0.5 Gyr, due to the fact that sub-mm 
output peaks half-way through spheroid formation. Indeed, depen- 
dent on assumed seed mass, the black-hole mass at the time of peak 
sub-mm emission is only 10 4 — 10 6 M©, offering an interesting al- 
ternative explanation (c.f. pure obscuration) for the failure to detect 
bright SCUBA sources at X-ray wavelengths (see Section B). 



4.2 The black-hole:spheroid mass relation 

Given the rather successful prediction of our fiducial model that a 
host spheroid with mass M Bp h — 10 12 M should house a black- 
hole of mass Mbh — 1.5 x 10 9 Mq, it is obviously of interest to 
undertake the same ca lculation for the range of spheroid masses 
originally modeled by Friaca & Tetievich (1998). The results of 
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Figure 6. Final black-hole mass plott ed against final spheroid m ass, for the 
range of galaxy masses explored by Friaca & Terlevich (|l 99q). The solid 



and open circles indicate the results assuming initial seed masses of 10 and 
1000 Mq respectively. The crosses indicate the 17 low-redshift galaxies 
which currently possess the most reliable estimate of black hole mass. The 
values of M a „h plotted for these objects are the revised values given by 
McLure & Dunlop ( 2002 ) from their reanalysis of galaxy bulge luminosity. 



Both this dataset and the model prediction are clearly consistent with a lin- 
ear relation of the form = k M sp /- l . The solid line indicates the best 
fitting relation to the model prediction, with a normalization k = 0.002. 



this calculation are presented in Figure |g, where predicted black- 
hole mass (for seed masses of 10 and 1000 Mq) is plotted against 
spheroid mass for five different galaxy masses. The best-fitting re- 
lation to these predicted data has a slope of 0.002, consistent with 
current estimates of the normalization of the Mhh '■ M xn h relation 

9- 



(Kormendy & Gebhardt 2001 



Merritt & Ferrarese 



2001 



The fact that this relation can be successfully predicted with- 
out any need to appeal to a causal connection between black-hole 
radiation and gas ejection is perhaps unexpected, but obviously in- 
teresting. It suggests that this relation may simply come about be- 
cause the amount of gas delivered to the central regions of a form- 
ing spheroid is roughly proportional to the total baryonic mass of 
the spheroid, with the rate and duration of gas delivery to the central 
region being limited by the star-formation process. 

4.3 AGN Activity 

Detailed observations of low-redshift AGN indicate that produc- 
tion of a QSO with My < —23 requires a black-hole mass 
Mbh > 10 8 Mq while production of a radio-loud source with 
P 2 . 7 GHz > W 25 WHz'V' 1 requires M bh > 10 9 M ( Dunlop 
et al. |2002j ; |Laor||1998| ; [VIcLurc & Duii1of|[2002| ). Although current 
constraints are poor, existing data on quasars at high redshift cer- 
tainly do not contradict the hypothesis that these mass thresholds 
are physically significant and still apply even at z ~ 2 (Kukula 



et al. ^2001 ). If we now apply these mass thresholds to the model 
predictions shown in Figure |5| we can see that the model predicts 
that quasar-level optical activity cannot be produced until ~ 0.5 
Gyr after commencement of star formation in the host, while pro- 
duction of a powerful radio source is delayed by ~ 1 Gyr. 
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However, the model predictions potentially extend beyond the 
basic prediction of such time delays. In particular, reference to Fig- 
ures q and ^| shows that, until the black hole passes the afore- 
mentioned mass threshold of ~ 10 8 Mq, the rate of gas delivery 
into the central regions of the galaxy by the cooling flow is well 
in excess of that required to sust a in Ed dington-limited black-hole 
growth. As discussed by Fabian ( 1999j ), in this situation it might 



10 



be reasonably assumed that all optical emission, and a substantial 
fraction of soft X-ray emission from the growing black-hole is ex- 
tinguished and that in its rapid growth phase the black hole may be 
only readily observable at hard X-ray wavelengths. 

The model thus makes a number of predictions about the ob- 
servability of high-redshift AGN at different wavelengths which 
can be compared with the results of current and future observations. 
These can be summarized as follows. 

(i) Assuming that, due to bias, massive spheroids such as that 
represented by our fiducial model are the first to collapse, the first 
AGN observable at optical wavelengths should be rather luminous 
QSOs, emerging ~ 0.5 Gyr after the onset of massive star forma- 
tion in the galaxy halo. This is because, by the time the central 
black hole ceases to be smothered by the central cooling flow and 
enters the regime of sub-Eddington accretion, it has already grown 
to a mass Mbh > 10 8 Mg. This prediction is at least consistent 
with the relativ e success of bri ght quasar surveys in locating high- 
redshift QSOs ( |FanetaL||200C| ). 

(ii) For a formation redshift z = 5, a powerful QSO should 
therefore emerge at z ~ 3.5, simply because Eddington-limited 
growth of the central black hole from a seed mass < 1000 Mq 
requires ~ 0.5 Gyr to produce a black hole of the required mass 
> 10 s Mq. For any formation redshift, luminous QSOs should 
not exist at redshifts significantly in exces s of z — 7. Th e current 
quasar redshift record stands at z — 6.28 ( ^an et al" 2001). 

(iii) For a formation redshift z — 5, production of a powerful ra- 
dio source is predicted to be delayed until z ~ 3, based once again 
on Eddington-limited growth and adoption of a minimum black- 
hole mass threshold of 10 9 Mq . For any formation redshift, pow- 
erful radio sources should not be found at redshifts significantly in 
excess of z = 5.5. The current radio-source redshift record stands 



at z = 5.19 (vanBreugeletal. 1999) 



(iv) Conversely, the predicted ~ 0.5-Gyr time lag between the 
onset of massive star-formation activity and the appearance of a 
central QSO in our model can be used to infer the first epoch of 
massive star-formation from the redshifts of the most distant QSOs. 
Given the current QSO redshift record of z = 6.28, our model 
would predict that massive star-formation commenced in its host 
galaxy at z ~ 10 — 12. 

(v) The very highest redshift AGN should only be visible in hard 
X-rays, as they rapidly approach a mass ~ 10 8 M via Eddington- 
limited accretion fueled by a cooling flow delivering gas at super- 
Eddington rates to the central regions of the forming galaxy. 



5 THE STARBURST-AGN CONNECTION 

Finally we focus on the issue of the connection between host- 
galaxy and black-hole evolution. In Figure ^ and Figure |§] we plot 
together the redshift dependence of AGN output as inferred from 
Figures ^ and ^| alongside the evolution of host sub-mm luminos- 
ity, again assuming a formation redshift Zf = 5. 

This combined diagram is instructive because it allows one to 
predict the properties of objects uncovered at different wavelengths 
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Figure 7. Predicted redshift evolution of the sub-mm output of a massive 
elliptical galaxy formed at Zf =5 compared with the predicted redshift 
evolution of the bolometric output of its central AGN, again shown for as- 
sumed black-hole seed masses of 10 and 1000 Mq respectively. The dashed 
line (left-hand axis) displays the logarithm of AGN activity, in terms of the 
fueling rate of the central black hole. The thick solid line marks the redshift 
at which the black-hole growth becomes sub-Eddington. The thin solid line 
(right-hand axis) shows the sub-mm output of the host, in terms of predicted 
observed 850 irm flux density. This plot demonstrates that for Zf = 5 and 
the currently-favored cosmological model, a massive forming spheroid can 
achieve peak AGN output at z ~ 3.5, while still being detectable in the 
sub-mm with Sgso — 10 m Jy, thus reproducing almost exactly the ob- 
served properties of the z = 3.8 radio galaxy 4C41.17. 



as a function of flux-density and redshift. We emphasize that al- 
though Figure [7] represents the predictions of the specific model 
considered here, most of the resulting observable consequences ap- 
ply to any quasi-monolithic model of elliptical formation in which 
star-formation activity peaks in a small fraction of a Gyr, while the 
rate of black-hole growth is dictated by Eddington-limited accre- 
tion starting from a modest seed mass. 

First, it can be seen that the predicted duty cycle of quasar ac- 
tivity at high redshift is considerably longer (~ 1 Gyr) than that 
which is generally inferred from (presumably retriggered) activity 
at low redshift, with excessive black-hole masses being avoided by 
the fact that most of the accretion during this primary, long-lived 
active phase is sub-Eddington. Consequently it is hard to escape 
the prediction that, for a reasonable range of formation redshifts, a 
substantial fraction of the present-day massive elliptical population 
should be active at the same time around z = 2 — 3. In fact, pro- 
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Figure 8. As Figure but with the AGN activity displayed on a linear 
scale, to better illustrate the rapid rise and demise of AGN activity as star 
formation in the host approaches completion. 



vided that one assumes at least half of optical QSOs are obscured 
by orientation effects, the comoving number density of quasars at 
2.5 is ~ 1 x 10~ 5 Mpc -3 , comparable to the present-day 



number density of ellipticals with L > 2 — 3L* ( frCochanek et al 
2001b. 



Second, it is clear that sub-mm observations of powerful AGN 
are unlikely to uncover the most luminous sub-mm sources, but 
should be capable of yielding detections for those objects observed 
in the overlap period when the black-hole has grown to large mass 
but sub-mm luminosity has not yet declined to undetectably low 
levels. Statistically, this is clearly most likely to be the case for the 
highest redshift AGN, and we would therefore predict that the sub- 
mm detectability of powerful AGN should be a strong function of 
redshift, approaching unity at z > 3. As already discussed this does 
indee^djH^pear^toJje^dTe^asejOT 

et al. |2001| ), and we would predict that the sub-mm detectability 
of optical and X-ray selected QSOs should be a similarly strong 
function of redshift. 

Third, it can be seen from Figure [7] that, most strikingly for the 
assumed seed mass of 1000 Mq, the formation redshift of 4C41.17 
inferred from its near-IR to sub-mm color in section 3, does indeed 
result in maximum AGN output at its observed redshift of z — 3.8 
as well as a sub-mm luminosity at this redshift in excellent accord 
with the observed value of 5*850 = 12 ± 0.9 mjy (Archibald et al 



Given that 4C41.17 appears to be one of the most massive objects 
in the universe, Figure W also indicates that the most luminous sub- 
mm sources are unlikely to exceed the sub-mm output of 4C41.17 
by more than a factor ~ 2, and thus sub-mm surveys are unlikely 
to uncover essentially any sources brighter than Ssso — 20 mJy. 

Finally, we would expect such extreme SCUBA sources, and 
indeed the majority of sub-mm sources uncovered by the current 
relatively shallow sub-mm surveys to be biased toward the epoch 
of peak sub-mm output. As illustrated in Figure ^ and Figure ^ at 
this stage (which occurs in our model only ~ 0.2 Gyr after the on- 
set of star-formation), the growing black hole is likely not only to 
be heavily obscured, but to have a black-hole mass (and AGN out- 
put) 10 3 — 10 5 times smaller than that of luminous quasars. This 
means that, regardless of their precise redshift distribution, the frac- 
tion of bright sub-mm detected sources which should also prove de- 
tectable at X-ray wavelengths with current facilities is small. This 
prediction is at least qualitatively consistent with current results 
which indicate that only ~ 10% of bright sub-mm sources show 
any detertable_s^gjts_ofAGN_ac^ 



ny \ 

et al. [200C| ; |Almaini et al j |2002| ; [Barger et al.||200l| ). We emphasize 
that while the growing black-hole at this stage may well be heavily 
obscured, the model presented here predicts that SCUBA sources 
are unlikely to be detectable with current Chandra surveys primar- 
ily because the black hole is simply not big enough, rather than 
by insi sting that the sources are Compton thick with < 1% scat- 
tering (Fabian 1999); even for unobscured AGN the deep Chandra 
observations presented by Almaini et al. (20021 are incapable of 
detecting a black-hole less massive than ~ 10 ' M Q at z ~ 3. The 
model therefore also predicts that AGN heating makes a negligible 
contribution to the luminosity of high-redshift SCUBA sources. 

In a future paper we intend to explore the effect of apply- 
ing our model to track the behavior of the batyons within a CDM 
merger tree leading to a massive (10 13 Mq) collapsed dark matter 
halo at z ~ 3. However, for now we simply note that the sim- 
ple model presented here may well offer a good approximation 
to the average behavior of the baryons within such a merger tree. 
This is because, for such highly-biased regions within the high- 
redshift universe, the merging of collapsed dark matter haloes of 
mass ~ 5 x 10 10 Mq into a final halo of mass ~ 10 13 Mq is com- 
pleted in a timescale < 1 Gyr, comparable to the timescale of the 
star-formation and black-hole formation processes described in this 
paper. 
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